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Abstract 

Implications stemming from the inclusion of non-perturbative, confinining effects, as 
contained in the Stochastic Vacuum Model of Dosch and Simonov, are considered in the 
context of a, hypothetical, quark-quark 'scattering process' in the Regge kinematical 
region. In a computation wherein the non-perturbative input enters as a correction to 
established perturbative results, a careful treatment of infrared divergencies is shown 
to imply the presence of an effective propagator associated with the existence of a linear 
term in the static potential. An equivalent statement is to say that the modified gluonic 
propagator receives contribution from a tachyonic ghost state, an occurence which is 
fully consistent with earlier such suggestions made in the context of low energy QCD 
phenomenology. 
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1. Introduction 

From the theoretical point of view, forward scattering at very high energies (Regge kine- 
matics) presents a situation where one can readily apply eikonal approximation techniques 
[1]. In this context, such processes provide grounds for exploring long distance properties of 
the underlying fundamental theory for the implicated interaction. For the particular case of 
QCD, long distance behavior constitutes a fundamental issue whose exploration is not only 
relevant to high energy processes but also to low energy phenomenology. 

In the present paper, we revisit the (idealized) problem of quark-quark 'scattering' in 
the Regge limit, which has been extensively studied within the framework of pQCD [1-5], 
with the aim to extend the aformentioned analyses in a direction which takes into account 
confining aspects of the theory. Specifically, we shall rely on the premises of the Stochastic 
Vacuum Model (SVM), pioneered by Dosch and Simonov [6] for the explicit purpose of 
accomodating the confinement property of QCD. A similar approach has been pursued by 
Nachtmann [7] by employing a different methodology from the one we shall adopt in this 
work. Different will also be the the scope of the present analysis. 

The construction of the SVM is motivated by the intention to incorporate established 
observations/results regarding the structure of the QCD vacuum [8] into a well defined 
theoretical framework. In particular, it summarizes all that is known and/or surmized about 
its properties through a set of three axioms, which are expressed in terms of field strength, 
as opposed to field potential, correlators. The underlying stochasticity assumption for the 
vacuum state facilitates the application of the cumulant expansion [9], which describes the 
factorization rules for higher order gluon field strength correlators in terms of two-point ones. 
One of the first results arrived at through the SVM is the deduction of the area law for the 
static Wilson loop, i.e. confinement. Specific applications of the SVM scheme, including 
comparisons with lattice results can be found, e.g., in Ref [10]. 

A concrete, as well as practical, way to apply the SVM scheme to specific situations has 
been suggested by Simonov [11]. The idea is to use the background gauge fixing method 
[12] and assign the background gauge fields with the task of becoming the agents of the non- 
perturbative dynamics. Specifically, one employs the gauge potential splitting A® = a® + B® 
with the a a ^ being associated with the usual perturbative field modes. The B°, on the other 
hand, enter as dynamical fields, assigned with the task of carrying the non-perturbative 
physics through field strength correlators which adhere to the cumulant expansion rules. 

Following our previous work of Ref [5] we find it convenient to employ the FFS-worldline 
casting of QCD, originally pioneered by Fock [13], Feynman [14] and Schwinger [15] and 
most recently revived in path integral versions, see Refs [16-18]. The reason for such a 
choice is that the eikonal approximation acquires a straightforward realization in this scheme, 
since, in the perturbative context at least, it can be readily implemented by restricting 
one's considerations to straight worldline paths. As it will turn out, the inclusion of input 
from the SVM will produce a kind of deformation of the eikonal paths, with low energy 
consequences, which represent, subleading, perturbative- nonperturbative interference effects. 
The conditions which justify the relevant computation will be made explicit in the text. 
Suffice it to say, at this point, that it was explicitly demonstrated in Ref [19] that the FFS- 
worldline formulation of QCD, in combination with the background gauge field splitting, is 
ideally suited for providing a framework within which one can directly and efficiently apply 
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the premises of the SVM. 

The main result of this paper is the following: Once (subleading) contributions incorpo- 
rating nonperturbative, as induced by the SVM, corrections to the perturbative expressions 
are taken into account, then a consistent analysis of the infrared issues entering the quark- 
quark high energy forward 'scattering' process reveals that the gluonic propagator exhibits a 
behavior which can be interpreted in terms of the presence of a 'tachyonic mass pole'. Argu- 
ments in favor of such an occurence, with important phenomenological as well as theoretical 
implications, have been promoted, from different perspectives, in several papers [20-22]. 

The organization of the paper is as follows. In the next section we present the basic 
formulas related to the 'amplitude' for the high energy 'quark-quark scattering process' in 
the forward direction and display their FFS-worldline form. Section 3 focuses its attention 
on infrared issues associated with the pertubative-nonperturbative interference effects under 
consideration in this study. The resulting expression for the amplitude is shown to be equiv- 
alent to the introduction of an effective propagator, which is associated with the presence 
of a linear term in the static potential. Section 4 extends the implications of the aformen- 
tioned result to issues related to renormalization: Modified running coupling constant and 
summation of leading logarithms via the Callan-Symanzyk equation. Finally, the technical 
manipulations leading to the main result of section 3 are displayed in the appendix. 

2. Preliminary considerations 

Consider an idealized quark-quark 'scattering' process in the Regge limit, defined by 
s/m 2 — > oo, s ^> t(— —q 2 = q± 2 ). The amplitude is given by 



T tlfjf (s/m 2 ,q 2 J\ 2 ) = J d 2 be-^ b E tt , ljf (s/m 2 ,l/b 2 \ 2 ), 



(1) 



where A is an 'infrared' scale, b is the impact distance in the transverse plane to the direction 
of the colliding quarks (assumed to be travelling along the rr 3 axis) while Eu'jj/, which 
incorporates the dynamics of the process, is specified, in Euclidean space-time, by [5] 



En',? = ( Pexp 



ig 



drvi ■ A{v\t) 



Pexp 



ig 



drv 2 ■ A(v 2 t + b) 



(2) 
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The vi and v 2 are constant four velocities characterizing the respective eikonal lines of 
the quarks participating in the 'scattering process' and 'conn' stands for 'connected'. In 
Minkowski space (and in light cone coordinates) one has v\ ~ (vf, 0, 0j_), v 2 — (0, v 2 ,0±), 
b ~ (0, 0, b±) with Vi ~ v 2 — -^^fsjm. 

Employing the gauge field splitting = + B a , the expectation values with respect 
to field configurations acquire the form < • • • >a=< • • • >a,B- Expanding in terms of powers 
of the perturbative field components, one obtains 




/oo r /'oo 

drv\- B{v\t) Pexp ig \ dr v 2 ■ B (v 2 t + b) 
-oo J ii> L J — oo 



33' I B 



ds\ ( P exp 



ig 



/ drvi ■ B{v\t) 

J si 



P exp 



ik 



ig 



/Sl 
drvi ■ B{vit) 
-oo 
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xP exp 



roo r rS2 

ig I dr t> 2 • B(v2T + b) Pexp ig / dr v 2 ■ B(v2T + b) 

JS2 J j m L J—oo 



P exp 


ig f 


- jm 


J — 



xt a kl t b mn v 2 ,v 1 ^G b ; u (v 2 s 2 + b,v 1 s 1 )) + 0(g 4 <a 4 >). (3) 



B 



Our objective in this paper is to study the behavior of the amplitude as |6| — > 0. Accord- 
ingly, contributions attributed exclusively to the non-perturbative, background terms will be 
ignored given that, by definition, they are finite in this limit. This narrows the expression 
of computational interest to the following one 



-g 2 N 2 _ ^ ii'tfj' J ds zj ds 1 V2 lx v lu (—Tr A iG lxu (v 2 S2 + b,v 1 s 1 )j 
+0(g 4 < a 4 >) + 0(b 2 ) (4) 



The propagator %G h ^ v (v 2 s 2 + b, ViSi) =< c^(t> 2 S2 + ^) a l{ v i s i) > acquires the following 
worldline expression 



poo r i nT o. / rT pT \ ba 

iG^(v 2 s 2 + b,v 1 s 1 )= j o dT J Vx(t)e~* Jo dti (t) P (e 9 ^ dtJ - F ^h dtiB j . (5) 



x(0) — vi s-^ 

x(T)=v 2 s 2 +b 



In the above formula (J-F)^ v = J^F a p, with jff the generators for the spin-1 representation 
of the Lorentz group 1 . It should also be noted that we have employed the notation B b ^ = 

B;(t c G ) ab = ~B c J ahc . 

To close this section let us briefly comment on gauge symmetry related issues. Given the 
'idealized' process under consideration our handling of gauge invariance will be to let the 
two worldlines of the 'colliding' quarks to extend to infinity in both directions and impose 
the boudary conditions — > as tEuci — ► ±oo. In this way, the overall worldline 
configuration introduces a Wilson loop in the path integrals, given that the end-points of the 
two trajectories can now be joined at +oo and at — oo. These, of course, do not correspond 
to boundary conditions for a scattering process per se, however our only objective in this 
work is to extract long distance implications based on the exchanges taking place in the 
immediate vicinity of the points of closest approach between the two worldlines. The study 
of realistic situations involving the scattering of physically observable particle entities in the 
Regge limit, using the presently adopted methodology, is under current consideration and 
the relevant analysis will be presented in the near future. Finally, concerning the issue of 
gauge fixing for the B field sector, our choice is prompted by the intention to rely on field 
strength correlators for describing nonperturbative dynamics [11]. It, accordingly, becomes 
convenient to work in the Fock-Schwinger (F-S) gauge [23, 24]. The latter is specified by 



rx pi 

a (x) = - du l/ (d ll u p )F°: u (u) = -(x - x ) u / daaF^xo + a(x - x )). (6) 
Jxq Jo 



1 Its incorporation into the worldline form of the propagator serves to signify the spin of the accomodated 
modes. 
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Of course, the arbitrary point Xq should not enter any gauge invariant expression. 



3. Infrared issues associated with the propagation of gluonic modes in a 
confining environment 

Consider the quantity defined by 

1 N c 
= N 2 _ ^ v^iTrAG^if^B = j^—jv 2fl v lv 



J™ dT J Vx(t)e-^o dt±2 ^/^Tr A Pexp(g£dtJ-F + igJ^ 




X / d r i 

JO 

x(0)=0 

x(T)=l 

where we have introduced 1^ = v 2 ^s 2 — v^Si — b^. It describes the propagation of the 
perturbative gluon modes in the presence of the background gauge field modes B a and, in 
this sense, it is expected to incorporate confinement effects associated with the SVM. 

Generally speaking, one would expect that in a study of a physically relevant process 
with the full (and proper) inclusion of non-perturbative effects, no need would arise for the 
introduction of an infrared cutoff to regulate the various expressions entering the computation 
at long distances. An infrared scale should, in other words, naturallly arise suppressing 
contributions from the very large distances (|/| — > oo). Given that the object of the present 
study is to investigate perturbative/non-perturbative interference, long distance effects in 
the (nonphysical) process of quark-quark high energy 'collision' in the forward direction, it 
becomes necessary to regulate infrared divergences associated with the upper limit of the 
T-integral. Our choice of introducing the infrared cutoff is via the relacement J °°(- ■ -)dT — > 
Jq 00 dTe~ TX (■•■)• O n a simple dimensional basis and given the length scales entering the 
problem, one could make the association A oc a\l\. 

Upon expanding the exponential one obtains 



x 



J Vx(t)e-*fo dt * 2 M ^ u l_ TrF{g±{t2) . Bixit^gxih) ■ B{x{t l ))) 1 



x(0)=0 

x(T)=l 



+ ^Tr F (gF, p {x{t 2 )gF; v {x{h))B} + 0((g^) B ) (8) 
Observe, now, that in the F-S gauge, the following relation holds 

Tr F {gB lta (x(t 2 ))gB ltl (x(t 1 ))) B = (x2-x ) va (x 1 -x ) vl J da 2 a 2 da x a x 

xTr F (gF^ 2l/2 (x + a 2 {x 2 - x^gF^^x^ + a 1 (x 1 - x ))) B , (9) 

which brings into play the field strength correlator. 
Setting Ui = x + aiix(t), i — 1, 2 one writes 

1 1 

-(9K 2 u 2 (u 2 )gF^ Ui ( Ul )) B = —Tr F (<P(xo,u 2 )gF^ U2 (u 2 ) 



2N C w"*^"* L " D n c 

X(f)(u 2 ,x )(f)(xo,u 1 )gF fllUl (u 1 )(f)(u 1 ,xo)} B = A { ^ U2 (u 2 - u x ) , (10) 
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where 0(xo,Wj) = Pexp (ig dv ■ B(vfj and is unity in the F-S gauge. Its insertion serves 
to underline the gauge invariance of the field strength correlator. 

With the above in place and upon making the redefinition U — > TU, i = 1,2, one deter- 
mines 

1(1) = ^ - j^v, ■ v 2 J o da 2 a 2 J o da iai J q dt 2 J q dt, 6(t 2 - h) 



x 



x(0)=0 

x(T)=l 



+ ^ 2 (^)^ 2 (t 2 )i Ml (ti)^ 1 x ( t 1 )A^ 2iAti ,Ja 2a ;(t 2 ) - a lX (ti)] \ + 0((g^) B ). (11) 



J^2 

On a kinematic basis, the correlator can be represented as follows [6,10] 

A M 2 k^i(^ - zi) = ((W'Wi - W^i) £) ( 2;2 ) 



+ 



1 9 



2 1 ^ r 2 1 



2 dZf J-1 i ■ 



One notices [10,11] that the first term enters as a distinct feature of the non-abelian nature 
of the gauge symmetry (it is not present, e.g., in QED). According to the premises of the 
SVM it is associated with the (QCD) string tension a by [10,11] 



Jo 



00 1 r 2 

dz 2 D(z 2 ) = - / d 2 zD(z 2 ) = -a. (13) 

7T J 7T 



Now, the central objective the present analysis is to determine first order contributions to 
the amplitude coming, via the SVM, from the non-perurbative/confining sector of QCD. 
The corresponding, lowest order correction is expected, on dimensional grounds, to be of the 
string tension a. This imlplies, as already pointed out by Simonov [25], that we shall set 
aside the Di term entering the kinematical analysis of the correlator, according to Eq (12), 
which cannot furnish contributing terms of dimension [m] 2 . In the appendix the following 
expression for the quantity 1(1) is established 



4tt2 W 



l + aa\l\ 2 a\n[C^) + 0(a 2 \l\ 4 ) 



(14) 



where a = | jV ^ 1 (l — k) with the constant parameter estimated to be k ~ 0.5. As pointed 
out by Simonov [21], the first -and most important- term contributing to a comes from the 
paramegnetic, attractive interaction between the spin of the gluons with the non-perturbative 
background field, cf. Eq (5). It should also be noted that the constant C entering the 
argument of the logarithm is connected with the choice of parametrization of D(z 2 ) see, 
eg., Ref [7]. In the context of a corresponding result having to do with an amplitude for 
a physically relevant, hence protected from infrared divergencies, process, then any such 
parameter would disappear. 

Suppose the following problem is now posed: Given the above results, which pertain to 
the gluon propagation in a confining environment, look for an equivalent, effective particle- 
like mode propagation, summarizing their full content. In this spirit, we shall proceed to 
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assess the possibility that the gist of all we have done up to here can be reproduced via the 
introduction of an effective propagator. Following Ref [20] we make the substitution 



1 1 E_ 

p ^ p + F 



(15) 



whose additional term signifies the presence of a linear term in the static potential. Then, 
one would obtain 



1(1) =v l -v 2 J 



j£k_ 



Sk-l 



(16) 



The integral is infrared divergent and should require the introduction of a corresponding 
cutoff. Alternatively, one could restrict the validity of the replacement, according to Eq 
(15), to the region k 2 > /i 2 . Then, one would determine 



7(0 



Vi ■ v 2 

k 2 >y? 
Vi ■ V 2 1 



d A k 

w 



ik-l 



JL_ fj? 

k 2 + k* 



1 + 



/x 2 |/| 



-ln- 



+ 0(//|/| 4 ) 



\l\ 2 (2vr) 4 [ 4 
Comparing the above result with that of Eq (14) one deduces 

3 N 2 



for/x 2 |/| 2 < 1. 



ttN 2 -1 



(1 - k)<t ~ 2.15a ~ 0.4 GeV 2 , 



(17) 



(18) 



in full accord with the phenomenologically determined estimate for the tachyonic 'pole' [20- 
22]. One also observes that A ~ \l\cr, as per our original supposition. 

Returning to the original, full expression, which provides the full dynamical input for the 
amlitude, we write 



Eaijji 



47T 2 I/| 



/-Too r- too 
ds 2 / ds x 
-oo J — oo 

r+OO r+OO I 

-g Ct a jj ,v 1 -v 2 J ^ ds 2 J ^ dsi^i 



■ V 2 



-oo 

"+00 



+oo 



l\ 2 



1 + 



n 2 \i\ 



■In 



/+oo f-\-oo P 

ds 2 \ dsi j 
-oo J ~oo J 



d A k 



'ik-l 



k 2 



+ 











4 




e/i 2 


i\ 2 


/ i 2 ^ 









(19) 



Concerning the logarithmic factors entering the above result, it is useful to remark that 
the various constants appearing in the arguments are not of any particular importance -at 
least in the approximation we have been working- given that they would disappear with the 
appropriate choice for the infrared cutoff. More importantly, thinking in terms of the signifi- 
cance of the above results if they became part of an amplitude corresponding to a physically 
consistent, hence protected from infrared divergencies, process, then any dependence from 
these constants should be absent. 

Going over to Minkowski space, the previous relation assumes the form 



Eu'jji ~ -g t^tjjdvi ■ v 2 / ds 2 / 

J —CO J— t 



dsi 



k 2 >fi 2 



d A k 



-ik-l 



k 2 + k 4 , 



(20) 
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Upon observing that 

/+oo r+oo 
ds 2 / d Sl e- ik - vlsl - ik - V2S2 = (2n) 2 v 1 -v 2 S(k-v 1 )S{k-v 2 ) = (2tt) 2 coth 7 5(A; + )5(A;_), 
-oo J — oo 

(21) 

where 7 is determined by 

cosh7 = 1 — n — r = =>- 7 ~ mis m ) =>- coth7 ~ 1. (22) 

|vi||f 2 | 2m 2 

It follows 

- "f^Wi'* coth 7 /(&V) (23) 



where 



/(&V) = ± / *k ± e*-+ (± + . (24) 



7T 

One immediately notices that if fi 2 b 2 <C 1, then 

/(/i2&2) * ln (^) ' (25) 

which recovers the known perturbative result -with an infrared cutoff given by A 2 = 
As b grows, while remaining in the region jj?b 2 < 1, one finds 

In turn, this gives 

E u , jf c -f^,^coth 7 {in (-*L) (l " ^) + • (27) 



4. Summation of large logarithms 

The presence of terms ~ g 2 lnl/b 2 fi 2 , entering through the function f(b 2 fi 2 ), imposes the 
need of their summation in the perturbative series. In the absence of the background field, 
i.e. in the framework of pQCD, it is well known that such a summation can be accomplished 
by employing the renormalization group strategies, which, for the quark-quark 'scattering' 
process under consideration, can be justified on the basis that b^ 1 plays the role of an 
ultraviolet cutoff. As Simonov has shown [11] the presence of the background field does not 
alter the, relevant for the summation, Callan-Symanzyk (C-S) equation. The physical basis 
on which this is so can be articulated by the following two arguments: 

1. Contributions from the non-perturbative sector do not introduce additional divergen- 
cies, given that they are finite at short distances. 
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2. Dimension carrying quantities arising from the non-perturbative sector (correlators) 
are structured in terms of combinations of renormalization group invariant quantities 
gB, i.e. they behave as external momenta, as opposed to masses which are subject to 
renormalization. 

Consequently, the called for renormalization group evolution follows the footsteps of the 
procedure employed in the purely perturbative analysis of the same 'process' [5]. In this 
connection it is recalled [2,26] that the Wilson contour configuration associated with E iV ^i 
mixes with one corresponding to a pair of closed loops resulting from an alternative way of 
identifying the points at infinity [5]. It is associated with 

Eij'ji' = - < ^c F ('ycoth'y-l)5 ij >5 j i>f(b 2 fi 2 ) 

+— c F (7coth7 - 1 -i7rcoth7)t",t",icoth7/(&V) + 0(al). (28) 

7T J J 

Accordingly, and upon introducing, in shorthand notation, W\ = 5u>5jj> + Ea>jj> and W 2 = 
Sij'Sji' + Eifji', the C-S equation assumes the form 



W a = -T ab W b , a, 6 =1,2, (29) 



with M playing the role of the uv cutoff whose running takes place between some lower 
scale (at which corresponding initial conditons are set) and an upper scale set by 1/6. The 
anomalous dimension matrix T ab , computed in the context of perturbation theory, see Refs 
[2,3,5], reads 

(r \ = ^1 ( ~f coth ^ incoth'j \ . 2 , . . 

1 ab) 7i { -icothi + 1 + mcoth-i N{coth-f - 1) - f coth-f ) + { sh 1 ' 

What does change, with respect to the perturbative analysis, on account of the the presence 
of non-perturbative, background contributions is the dependence of the W a on the -B-field 
correlators, i.e. one has W a = ^[{A^}, M,g}. Given that the computation has taken into 
account only the two-point correlator, the extra dependence of the W a will involve the string 
tension, cf. Eq (13). The inclusion of this additional dimensional parameter will have its 
effects on the running coupling constant. 

With this in mind, let us recast Eq (29) in integral form: 



W a [<T,M 2 ,g B (M 2 )] = {Pexp 



M 2 dM 

^T(g B (M)) 

Mi M 



> W b [a,M 1 ,g B (M 1 )\ (31) 

ab 



with the path ordering becoming necessary because the anomalous dimension matrices do 
not commute with each other. Concerning the integration limits a consistent choice, given 
the premises of the present calculation, is to take M 2 = 1/6 and set M 1 = l/b with b\a < 1. 
It is observed that the non-perturbative input enters the W a not only through their explicit 
dependence on the string constant, but also -which is the most impotant- through a running 
coupling constant g B , which obeys the equation 

m -^9b(M)=P(9b(M)). (32) 
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The solution of the latter calls for initial conditions which are influenced by the presence of 
the non-perturbative background and specifically by a. Such matters have been studied by 
Simonov in [25]. 

Turning our attention to the 'deformation' (to the one-loop order) of the running coupling 
constant, on account of its additional dependence on the background field, we proceed as 
follows. Knowing the perturbative result to order as, we go to Eq (31) and present its 
solution in the form 



W\[a, 1/6, Q!b(1/6)] = SaiSjji—i coth^f 
It follows that 



a s (l/bl)f(by) 1^ ^a s 



r 



tl>t%+0(a s ). (33) 



r LI o H T 

J -Mr) = a s (l/b 2 )f(b^ 2 ) - a s (l/bl)f(btn 2 ) + 0(a 2 s ). 



Upon comparing with Eq (27) one obtains 



u 

a B (r) = a s (r) il + — 



(34) 




+ 0(a 2 ). 



(35) 



It should be noted that the validity of the above results holds for r/fj 2 > 1 and as{r) 



47T 



0- i n ( T M2) < 1- An indicative estimate, on the basis of Eq (35), is that if as — 0.5, then 
as — 0.5(1 + 0.05). Following Refs [3,5], one surmises that the 'amplitude' A for the 'process' 
under consideration behaves as 



cxp 



■— -In — / — a B {r) + 0(a 2 s ) 



oc exp 



a s 



2tt 



^N c ln[ — )f(b 2 ^ 



(36) 



from which one 'reads' a reggeized behavior for the gluon. The difference from the usual, 
purely perturbative result is that the function f(b 2 fi 2 ) is now connected with the modified 
propagator, as per Eq (24). 

In conclusion, we have demostrated that the non-perturbative input, through the SVM, 
to the analysis of a hypothetical quark-quark 'scattering' process in the Regge kinematical 
region, produces a result which, in a phenomenological context, has been argued to be 
extremely attractive in reproducing low energy hadron phenomenology. In a sense, this 
investigation could be considered as a special example, which justifies Simonov's more general 
argumentation [22] according to which the perturbative-nonperturbative interference in static 
QCD interactions at small distances imply the presence of a linear term in the potential. 
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Appendix 

Given the set of the defining, worldline formulas given by Eqs (6)- (15) in section 3, we 
proceed to derive derive Eq (16). In the course of the derivation the various quantities and 
parameters appearing in the last part of the section are specified. 
Writing 

D(z 2 ) = / dpD(p)e- pz (A.l) 
Jo 

one determines 

1(1) = — 2 w Nc2 C _^ -V2j o da 2 a 2 J Q da iai J o dt 2 J q dt^-t,) 

x dpD(p)[48Q(p;t 2 ,t 1 )-R(p;t 2 ,t 1 ,a 2 a 1 )] + 0(<g 4 F 4 )> B ), (A.2) 
Jo 



where we have made the change ti — > Tti, i — 1,2 and have introduced the quantities 
Q(p-t 2 M)= f dTe ~ TX2 fj^ e ~^ J Vx^e-^fo^e^M-*^ (A.3) 



(0)=0 

x(l)=l 



and 



/ 7T \ f°° 2 f d Q i 

R(p;t 2 ,t 1 ,a 2 ,a 1 )= I -J dT e~ TX J j—^e~^ (5^5^, - 5^5^) 

x J Vx(t)e--^o dt±2 ^x^(t 2 )x V2 {t 2 )x^ (A.4) 

x(T)=l 

with \l\, as defined in the text. 

The above path integrals can be executed by employing standard techniques, given that 
'particle' action functionals are quadradic (plus a linear term) [16-18]. Ignoring terms giving 
contributions 0(b 2 ) and using condensed notation from hereon, one determines 

^khr^fw?-*-™-^™ (A - 5) 

and 

fl=~f dTe ~ T " 2 l -Pv^'^^ + T ^ + T Vc 2 + +0(l 2 q% (A.6) 
16 p z Jo J (2vr) 4 

where the following, one dimensional particle propagator-type quantities have been intro- 
duced 

Ai 2 = A(t u t 2 ) = fi(l - t 2 )9(t 2 - h) +t 2 (l- t^Oih - t 2 ), (A.7) 



G 12 = G{t 2 ,h) = A(t 2 ,t 2 ) + A(ti,ti) - 2A(t 1 ,t 2 ) = \t 2 - txKl - |* 2 -^1) (A.i 
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and 

K 12 = K(t 2 , h) = a 2 2 A(t 2 , t 2 ) + a\A{t u h) - 2a 1 a 2 A(t 1 ,t 2 ). (A.9) 
The coefficients entering Eq(A.6) are given by the expressions 

c = -72A 12 9 1 A 12 9 2 A 12 (A. 10) 

ci = 48Ai 2 (9 2 Ai 2 (9i^i 2 + 24ai<9iA 12 (9 2 Ai 2 (aiAi 1 - a 2 Ai 2 ) 

+24a 2 (9iAi 2 i9 2 Ai 2 (a 2 A 22 - «iAi 2 ) - 12(alA 12 diA u d 2 A 12 + a 2 Ai 2 <9 2 A 22 <9iAi 2 
-aia 2 Ai 2 9iAn(9 2 A 22 - aia 2 Ai 2 <9iAi 2 <9 2 Ai 2 ) 
-12(a 2 A 22 - a 1 A l2 )(a 1 d 1 A ll d 2 A 12 - a 2 <9iAi 2 <9 2 Ai 2 ) 

-12(aiAii - a 2 Ai 2 )(a 2 (9 2 A 22 (9iAi 2 - ai<9iAi 2 <9 2 Ai 2 ) (A.ll) 



and 



c 2 = 24(a 2 A 22 - « 1 A 12 )(a 1 A 11 - « 2 A 12 )9 1 ^i 2 9 2 A 12 . (A.12) 
Given the above, the "paramagnetic" contribution to Eq (A. 2) becomes 



2N r°° ~ r 1 r 1 

12 at2 c v 1 -v 2 dpD(p) dt 2 dhOih-t^Qip-t^ti 
Nc A — 1 jo jo jo 

12 2Nl_v 11 v 1 rdp b{ ) 1„{ A „-^P \ , ffil2 

1 4vr 2 Jo v 



ln(4e-^^) +0(A7p) 



(A.13) 



and since 



Eq (A.9) gives 



16JV C -1 4tt 2 Jo p 

r^bip) = r 'dz 2 b{z 2 ) = - r d 2 zb{z 2 ) = ~ a , (a.u) 

JO V JO 7T JO 7T 



7 ' = jV?-l 4ir* ffln l g ^j- (A ' 15) 



This furnishes the correction term from the background gauge field contributions entering Eq 
(14) in the text. The constant C entering the above result depends on the parametrization 
of D(z 2 ). Following the one of Nachtman [7], one determines C = 39.65. The numerical 
computation of the factor k, based on the expressions for the Cj, as given by (A.ll), produces 
the value k ~ 0.5. 
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